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We analyze the two-dimensional spin-fermion model in the strong coupling regime relevant to 
underdoped cuprates. We recall the set of general sumrules that relate moments of spectral density 
and the imaginary part of fermion self-energy with static correlation functions. We show that two- 
pole approximation of projection method satisfies the sumrules for first four moments of spectral 
density and gives an exact upper bound for quasiparticle energy near the band bottom. We prove 
that non-crossing approximation that is often made in perturbative consideration of the model 
violates the sumrule for third moment of spectral density. This leads to wrong position of lowest 
quasiparticle band. On the other hand, the projection method is inadequate in weak coupling limit 
because of approximate treatment of kinetic energy term. We propose a generalization of projection 
method that overcomes this default and give the fermion self-energy that correctly behaves both in 
weak and strong coupling limits. 

Short title: Generalized projection method 
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I. INTRODUCTION 

Over the past few years it became increasingly clear that anomalous normal and superconducting state properties 
of high- Tc cuprates are governed by their close proximity to the transition into antiferromagnetic Mott-Hubbard 
insulator state. The transition do occurs for strongly underdoped systems. There are many indications that even 
overdoped cuprates are never far from being antiferromagnetic Such kinds of incipient antiferromagnet at low 

temperatures have both fermionic and spin excitations that strongly interact. For description of such interaction the 
two-dimensional spin-fermion models are used. They have a general form 

H — Hkin + Hs + Hint- (1) 

Here the kinetic energy term 

p r 

describes the bare fermion propagation. We take for definiteness the simplest version of nearest neighbor hopping 
Cp = —2t{cospx + cospy). 

The spin subsystem may be described by microscopic Heisenberg Hamiltonian 
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or quite generally by a phenomenological form of dynamic susceptibility |^ we are interested by only 

fermion spectrum. In the present work we suppose that spin subsystem is in the paramagnetic state {S") = with 
strong antiferromagnetic correlations Cg — {S"_^_gS") < ,the correlation function Cq — N^^ '^^exp{iqr) {Sp_^_^Sp') 
is strongly peaked at wave vectors in the vicinity of Q = (tt, tt). 
The third term in its simplest form describes Kondo interaction 

r ^ pq r 

In the above formulae the sums run over the sites r of a square lattice and over the nearest neighbours with the 
lattice spacing | g |= 1. For short we miss the spin index for creation a^^ and annihilation a^^ operators of the 
Fermi particles (we shall name them as electrons) and in the Hamiltonian of Kondo-interaction Hint we take the 
notation Sr = S!^a°'] summation over repeated indexes is understood everywhere; cr" are the Pauli matrices. (...) 
means thermodynamic average over grand canonical ensemble. 

In order to achieve the more realistic description, the interaction term may be generalized by addition of nearest 
neighbour coupling a\.Srar+g + h.c. , etc. It is essential that in any case Hint remains local in real space and does 
not couple directly the excitations on sites separated by large distance. 

The usual approaches to the Hamiltonian (|l|) exploit its apparent similarity to that of polaron problem. Various 
perturbative methods have the advantage of exact treatment of the one-particle part, Hkin- It is widely believed that 
non-crossing approximation is appropriate for spin-fermion systems even in strong-coupling regime J ^ t . 

Below we shall prove that in this regime the non-crossing approximation violates the sum rule for the third moment 
of spectral density, as a result, it gives the wrong position of lowest 'singlet' band. An alternative is the Mori-Zwanzig 
projection technique due to the local nature of Hint it suffices small number of basic operators for appropriate 
account of local correlations. The obvious disadvantage of the technique is that the kinetic energy is treated in 
approximate way, as a consequence, it fails to correctly reproduce the weak-coupling regime J <ti t. In the present 
paper we exactly take into account the kinetic energy term and use the projection method for remaining terms. In 
the end we obtain the fermion self-energy that correctly behaves in the weak-coupling limit and gives the right lowest 
band position in strong coupling limit. 



II. SUM RULES FOR GREEN'S FUNCTION AND SELF-ENERGY 



The quantities we calculate are the retarded fermion Green's function (GF) 

fOO 

GxYiio) = (X|yt) = _^ die'"(*-*')({x(i),yt(t')}) . (3) 



and the spectral density 



Axy{^) = --Ini [GxY (w + «0)] 

TT 



Here and below {...,...},[...,...] stand for anticommutator and commutator respectively. For diagonal GF Y = 
X, the spectral density is positively definite AxxJ.^) > 0. Mori-Zwanzig projection method allows to represent 
Gxx{z),lTaz > in the continued fraction form 

b2 h2 h2 
Gxxiz) = ^ ^-^ (4) 

z — ao— z — ai— z — an— 

where 

The coefficients b„,a„,n > are related with the spectral density Axx{'-^) via the set of orthogonal polynomials 
Pn{Lo), satisfying the recurrence [p|,p^: 

P_i(c.) = 0, Po(c.) = l, 

Pn+l{uj) = {u;^ &n)Pn{^) - hlPn-l{uj), (6) 



bo = / Axxii^)duj, ^0 ^ -r2 I LoAxx{^)duj- (5) 
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and 



/-co P^+i{uj)Axxiuj)d^ 

Here we have used the nonnormahzcd form of the polynomials Pn{uj)Ps{i^')Axx{i^)duj — (^nsdlm^i t>m)^- 
On the other hand, from the equation of motion 

uj{X\Y^) = {{X,Y^) + {XC\Y^), XC=\x,h] (9) 

follows the sum rule 

F{uj)AxY{^)duj = {{XF{£),Y^) (10) 

for arbitrary function F{C). In particular, it establishes the relations of coefficients a„,b„ with static correlation 
functions 

hl = {{X,X^), ao^ ^lg'^y , (11) 

({X(£-ao)^Xt}) _ ({X/:(/:-ao)^Xt}) 

({X,Xt}) ' ({X(/:-ao)2,Xt}) ^^'^ 

Now introducing the self energy E(z) through the relation 

/ {{xc,xn) 



-j:{z)\Gxxiz) = {{X,X^) (13) 



and comparing (^ and ( p^ we see that S(z) it is the continued fraction similar to G{z) . Thus we can introduce the 
spectral density 

p{lu) = -Im[I](tj + iO+)]/tt 
and obtain for it sum rules that follow from ( p^ ) 

hl= p{Lu)duj, (14) 



1 /■ + °° 

ai = — 2 / ujp{Lj)duj. 



(15) 



As it follows from (|T^ and (||) the last equality relates the third moment of Axxi^) , the first moment of p{uj) and 
static correlation functions. For spin-fermion models in the limit of low doping, the spin-spin correlation functions are 
only involved in (^2|). Below we show that non-crossing approximation violates sumrule (15) and is obviously wrong 
in strong coupling limit J 3> i. 
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III. PROJECTION TECHNIQUE 



In practice, projection technique calculations are possible only for finite basis set and only finite number of continued 
fraction levels is possible to calculate in (^). 
In the frameworks of our model we have 

{u;-ep){ap\al) = l + J^2{bp\al), 6p = iV-^/^ ^ 6,e-'f^ (16) 



v/2 ^ ' 4 

Thus the 'bare' electron operator with one-site spin polaron operator br represent the natural basis set for appro- 
priate account of local correlations. It is important that this set is closed relative to Hint , i.e. 



Now the commutation relation 



J\/f2br, 



67-, Hint 



J V /2ar 



-t ^ ^ Sr+Rdr+g + J Sr+RSrO-r + Sr+Rdr, Hg 



(17) 



(18) 



is projected onto the basis set in order to decouple the equation of motion for higher order GF Gb{p,uj) = (6p|aJ) , 
In the following, we neglect the spin excitation energy I <^ t, J 



LoGh 



bp,H 



Ca 



— J ) Gh + JV /2Ga 



This gives both GF in two pole approximation 
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1 S |2 


1 T 1 


Km 4 





1, 



here f2„, ,i — a,b,n = S,T, (fig < fix) are eigenvalues and eigenvectors of the problem 

= 0, 



where matrix elements are 



ao — r^n bi 
bi ai — ilr, 



ao = Cp, bi = J\/j2: ai 



/2 



J. 



(19) 



(20) 



(21) 



(22) 



Together with the normalization bg = 1 the matrix elements correspond to continued fraction coefficients of g'^P that 

(2) 

coincide with the first two pairs of coefficients of exact GF Ga- It means that Ga automatically satisfies the sumrules 
(0), (ll)and dH). 
Near the band bottom Ga and Aaa should have the form 



Ga{p,uj) = 



Zgjp) 
LU — Eri 



Gi. 



(23) 



^aa(P, w) = Za{p)S{uj - Ep) + A,„c{P,i^)- 

Here Ep and Z{p) < 1 are quasiparticle energy and pole strength respectively. The incoherent part. Ainc is not zero 
for uj > LUmin > Ep. Now it is easy to show that fls represents an exact upper bound for Ep. Let us consider eigen 
operator for lowest 'singlet' band 
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The GF 

Gee = ^ + Gee.mc 

has the pole at the same energy as bare fermion GF Ga (in our model 7^ for all p). On the other hand, from 
and ( ^if) we have 



= ^e(p)^P + / w^e?,mc(p, = z^{p)E^, + [1 - Ze(p)] f)„c, (24) 

here 

^Hnc — noo 1 f y. J — "^min 

is the center of gravity of the incoherent part. As the pole strength is < Ze(p) < 1, we have 

Ep<ns< 17,„„ (25) 

i.e. for any p the exact energy Ep lies always lower then the energy given by two pole approximation. 
In the strong coupling limit that gives 

ns^-p + ep{^ + Cg), (26) 

and from (|2^) it follows that actual pole position is lower than Cls- It is not difhcult to calculate next continued 
fraction coefhcient 



in the approximate equality we took into account that main contribution to the sum over q comes from the vicinity of 
Q and N"^ Cq = f2- We see that b2 is of the order of kinetic energy t <^ J . It means that small polaron formed 
by our basic operators interact with the spin subsystem much weaker than the bare hole. So, the expected polaron 
energy renormalization from to Ep is of the order of / J . 



IV. EXACT TREATMENT OF THE KINETIC ENERGY 

From the above consideration it follows that the two-pole expression ( pO| ) is useful in strong coupling limit. In the 
opposite case J <C t it becomes not appropriate because it gives the self-energy in the one pole form 

I](2)(z) = b2/(z-ai) (28) 

and cannot describe the damping of quasiparticle. The reason is in the approximate treatment of the kinetic energy 
term in course of the projection of the equation for Gb ( |l^ ). Here we propose a generalization of projection technique 
that completely removes this default and makes it possible to exactly take into account the kinetic energy term. 
We express bp in the following form 

the equation of motion for every item gives 

{LU-ep+q)(Sqap+q\al\^{Sqap+q,V |at\. (30) 



Here and below we denote V = Hint + Hg. Now we project the higher-order operator in the right-hand side of (^ ) 
onto our basis operators Bi^p = Up, B2,p = bp 
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({ Sqap+q,V ,bI^Y/ = ^^^j= ^ (|| Sr^ar^,V ,B\^^Yieyi^{iqri-i{p + q)r2+iprz) 



— y({ 



+ »S'ri 



' } ) *5XP (''I - ?-2) + «P (J-S - ''2)) 



ftr, V 



,-B|^|\exp (iqR) = ^= V ifi.fl exp = -^K^^g, 

' J ' V jv ■^r V -iV 



(31) 



1 = 1,2 . We have used the local character of operator V that gives r^- Thus Eq.(|30|) may be rewritten as 



Now equation for bp is 



(32) 



(33) 



Explicit calculation gives (/ « 0) 

Kl^q ^ JCq, K2,q 

and we obtain the following form for the fermion self-energy 



J 



Co 



T.{p,uj) = — 'LJ^ — n _ 



C„ 



(34) 



J + hDp 

In the weak coupling limit J this expression coincides with the second order result of perturbation theory 

Epert (p, t^) = -'"^-Dp(w), (35) 



which corresponds to the projection ( |31| ) only on the first operator Op. In this limit a polaron of large radius is formed 
and bare fermion represents slightly damping quasiparticle. 

Now let us show that the self-energy (|3^) also give correct result in the strong coupling limit. The pole position is 
given by i?p — ep — S (p, /?p) =0 . For determination of lowest band position we may neglect the Cp+g oc t compared 
w oc — J and write Dp(io) « , then 



S (p, uj) 



JV2 



J + W/2' 



-C/n — 



The perturbation theory result is 



Spert (P,t^) ~ J^/2/'^, 



pert 



2>^ -J\lj>i^S- 
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The reason of the perturbation theory fail is the violation of the sum rule ( nsh . We have 



j: 

Ppert(p,t^) = — ^CqS{uJ - ep+q), 



N 



1 

^i.pert = :^J^^Cq J UjS{uJ - ep+q)duJ 



f2N 



h 



(36) 



Comparing (|3^) with the exact value given by Eq. (|2^) we see the absence of terms proportional to J that give 
completely wrong result for J ^ t. It is not difficult to prove that summation of the infinite series of non crossing 
diagrams for the self energy does not change the value of ai ^pert (^6|). Indeed, the non-crossing (self-consistent Born) 
approximation gives 



j- 



N 



(37) 



where ujq is the energy of spin excitations. We then have 



ai > 



^-^^Cg J UjAaa{p + q.i^-UJq)duj = j-^^Cq{epJ,q- 



■UJa 



and obtain the same result ( |3q) because Uq is negligible, at least in the vicinity of q = Q. 

Moreover, the self-energy ( pTl) leads to the absence of quasiparticles. We may write I]„_c(p, cj) « J^f2Ga {p + Q, w), 
then 



Ga{p,U}) = [w - ep - J^/2G'a(p+ Q,CJ)] 



Solution of the quadratic equation gives the expression for Green's function 



Ga{p,Uj) = 



[uj - ep){uj - Ep+q) - ^[uj - ep)2(w - ep+qY - AJ"^ f2{uj - ep)[uj - ep+g) 



2J2/2(c^ - ep) 

which has no poles. Analogous result was obtained previously in Ref. [Q. 



V. NUMERICAL RESULTS 



On Figure la,b we present the spectral density Aaaip,!-^ + iy), V — O.OSi that corresponds to three different repre- 
sentation of fermion self-energy: two pole approximation , Eq. (|2^), generalized projection method, Eq. (|3j), and 
second order perturbation theory, Eq. (35). We took the value J/t = 3 that is typical for underdoped cuprates [0-^. 
For spin-spin correlation function we took the expression 



Ca 



l 3\Cg\il-Jq) 

2ai(l + 79) ' 



which is provided by spherically symmetric theory for Heisenberg model on square lattice 1 11 ij] (7^ = (cos + 
cosqy)/2,Cg ~ — 0.35,q;i « 2.35). We calculate the function Dp{oj) by direct summation over n x n q— points 
in Brillouin zone (the results for n = 32 and n = 64 are almost indistinguishable). From Fiq.l we see that the 
lowest pole position obtained by generalized projection method satisfy the relation (Eq), in contrast to that given by 
perturbation theory. 

Figure 2 shows the spectral function in generalized projection method for various values of p along the diagonal 
of Brillouin zone. Quasiparticle pole exists throughout the whole Brillouin zone. Quasiparticle dispersion Ep has 
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only approximate symmetry relative to the boundary of antiferromagnetic Brillouin zone. Strong asymmetry of the 
spectral weight in the generalized projection method is the consequence of the sumrule (|]) - center of gravity Aaa{p, ^) 
should coincide with ag = ep. So, near the point p = (tt, tt) , where the quasiparticle peak is far from ep , its weight 
is small. 

We have shown above that the expression (^) gives reasonable solutions both in weak and strong coupling limits. 
Thus we may expect it to be valid in the intermediate regime J ^ t. Figure 3 demonstrate the fermion spectral 
density for J = \.bt . In this regime the quasiparticle solution exists only near the two band minima p — {tt,tt) and 
p = (0, 0) . In other points (for which we could say that quasiparticle pole lies within the band of bare fermions) the 
interaction mixes solutions with different p. For p = 0.2(7r, tt) and p = 0.8(7r, tt) we have resonant solutions near the 
bottom (emin = —4:t) and the top (cmax = 4i) of bare fermionic band respectively. Near p = (7r/2,7r/2) we have pure 
incoherent spectrum. 

VI. CONCLUSION 

We have considered spin- fermion model that is often used for description of strongly correlated systems. We have 
compared two popular approaches to the calculation of fermionic Green's function: the non-crossing approximation 
of perturbation theory and Mori-Zwantzig projection technique. We have shown that first one is valid only in the 
weak coupling regime and the second one only in the strong coupling regime. The non-crossing approximation gives 
wrong position of lowest quasiparticle band. The reason is the rough violation of the sum rule for the third moment of 
spectral density. The summation of the infinite series of non-crossing diagrams for the self-energy does not alter this 
result. We have proposed the generalized version of projection method that treat exactly the quadratic on fermion 
fields the kinetic energy term. The resulting expression for the self-energy coincides with that of perturbation theory 
in weak coupling limit and provides right quasiparticle pole position in strong coupling limit. We thus consider it as 
a good starting point to the investigation of intermediate coupling regime that is believed to be relevant for optimally 
doped cuprate compounds. 
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FIG. 1. The spectral density of the one particle Green's function with the self-energy, obtained by three different ways: two 
pole approximation (Eq), generalized projection method, (p4), and second order perturbation theory, (p5|) 
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FIG. 2. The spectral density in generalised projection method as a function of quasimomentum p in strong coupling regime 
J/t = 3. 

FIG. 3. The spectral density in intermediate paremeter regime. 
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